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Abstract 

Experimental results indicating the existence of the high magnetic field 
Wigner Crystal have been available for a number of years. While variational 
wavefunctions have demonstrated the instability of the Laughlin liquid to a 
Wigner Crystal at sufficiently small filling, calculations of the excitation gaps 

O 

have been hampered by the strong correlations. Recently a new Hamiltonian 
formulation of the fractional quantum Hall problem has been developed. In 
t — ' this work we extend the Hamiltonian approach to include states of nonuniform 

in ' 

density and use it to compute the excitation gaps of the Wigner Crystal 
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states. We find that the Wigner Crystal states near v = 1/5 are quantitatively 
well described as crystals of Composite Fermions with four vortices attached. 
Predictions for gaps and the shear modulus of the crystal are presented, and 
q \ found to be in reasonable agreement with experiments. 

> 

^ ■ I. INTRODUCTION AND PREVIEW 

.«?; 

The Fractional Quantum Hall (FQH) regime presents us with the quintessential problem 
of strong correlations. In a strong magnetic field B, the kinetic energy of the two-dimensional 
electron gas (2DEG) is quantized into Landau levels with energy (n+|)u; c , where u c = eB/m 
is the cyclotron frequency. Each of these Landau levels (LLs) has a huge degeneracy equal 
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to the number of quanta of flux penetrating the 2DEG. When the lowest Landau level 
(LLL) is partially full, it is seen that the kinetic energy is degenerate, and all the dynamics 
has to come from interactions. At certain special filling factors (recall that filling factor is 
v = 27m/ eB) the system reorganizes itself into new strongly correlated ground states |I[ 
with fractionally charged excitations 0. 

The past decade has seen the development and acceptance of the Composite Fermion 
(CF) concept as basic to the understanding of a variety of these electronic states ||. The 
CF is pictured as an electron bound to an even number / of quanta of statistical flux, which 
are opposed to the external field. At the mean field level, each CF sees both the external 
field and the statistical field due to the other particles, and therefore moves in an effective 
field B* = B — 27rln, where n is the density of electrons. The principal fractions v = p/2p + l 
are seen to be exactly those fillings when the number of particles is exactly enough to fill an 
integer number of LLs of the effective field. 

Thinking in terms of CFs greatly simplifies the description of different incompressible and 
compressible FQH states. CFs are believed to be the true quasiparticles in much the same 
way as Landau quasiparticles are for the normal Fermi liquid. Following Laughlin's seminal 
insight 0, the original CF theory was based on generating electronic wave functions [|J. 
Contemporaneously, field-theoretic approaches || were also developed to better understand 
the FQHE, and to compute response functions. Most of the field-theoretic approaches are 
based on the Chern-Simons (CS) transformation, a method of attaching flux to particles in 
two dimensions. Attaching an odd number of flux quanta to electrons transforms them into 
bosons and leads to the bosonic CS approach , while adding an even number leads to the 
fermionic CS theory |J. In the mean field approximation the fermionic CS theory recovers 
the picture of CFs in an effective magnetic field. Recently, based on the fermionic CS theory, 
a Hamiltonian approach was developed ]7[ to describe liquid states in the FQH regime. In 
this approach, the CF representation is reached from the bare electronic coordinates by 
a series of canonical transformations. The end product is the electron density operator 
reexpressed in the CF coordinates suitable for further calculations and/or approximations. 
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Physical quantities calculated in this approach seem to be in reasonable agreement with 
numerical results and experiments ||. 

The subject of this paper, however, is not the liquid FQH states, but the insulating 
states that have been detected experimentally at very low filling fractions ||. A natural 
candidate to exhibit such insulating behavior is the electronic Wigner Crystal (WC). The 
simplest description of this state is the Hartree-Fock (HF) wave function |10| 

i 

where A is the antisymmetrization operator, and 0^ is a single-particle wave function that 
is localized at (lattice site) and belongs to the LLL. It is given by 

where Iq = (eB)" 1 ^ 2 is the magnetic length. The wave function ([]]) has been improved 
by adding a Jastrow correlation factor fl2|| , and the energy of the resulting state has been 



shown to become lower than that of the liquid state at about the experimentally right filling 
fraction {y « \) 0-0. Thus, a very strong magnetic field favors crystalline order by 



localizing the electrons. 

However, not all the experimental evidence supports the simple electronic WC picture. 



In particular, transport experiments 113| , |14| suggest that the activation gap in this system 



is two orders of magnitude smaller than the theoretical estimate as calculated using the 



Hartree-Fock approximation [|T^,|T^]. Moreover, close to the Laughlin fractions 
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a dip in the longitudinal resistivity p xx is observed |L7j, resembling the behavior of the 
correlated liquid state. The measurements of the Hall resistivity p xy are surprising as well 
|r?Hl9[. The electronic WC is known to have a vanishing Hall conductance a xy = 0, which 
implies a vanishing Hall resistance p xy = 0. On the contrary, experiments see Hall insulating 
behavior, that is, p xy ~ These problems led Yi and Fertig |?(J to consider crystalline 
states with Laughlin- Jastrow correlations. The idea is to construct a trial wave function 
with correlation factors that keep electrons apart. Each electron is combined with I vortices 
to obtain the trial wave function 



=^n(*-*i)'n^( r <)- ( 3 ) 

*/j * 

The Coulomb energy for this wave function was then computed using Monte Carlo simu- 
lation. Yi and Fertig showed that the ground state energy of the correlated WC is lower 
than that of the usual WC at experimentally relevant filling fractions [^]. Moreover by 
introducing Laughlin-Jastrow correlations between the interstitials and the lattice electrons 
the experimentally observed p xy (Hall insulating behavior) ]21 , 22| can be explained. Un- 
fortunately, the method becomes too computationally demanding to allow one to calculate 
other quantities of interest, such as the excitation spectrum. 

Since the Laughlin-Jastrow correlations are precisely the ones that convert electrons into 
CFs [|J we are led to consider a crystal of Composite Fermions. The main advantage of 
the Hamiltonian approach is that one can easily compute the excitation gap along with the 
ground-state energy. 

Our main result is that the CF theory with I = 4 zeros attached to each particle gives 
the best description of the experimental phenomenology. Figure 1 shows the results of our 
calculations of the excitation gap E g as a function of the filling factor around v — 1/5 in the 
triangular lattices of CFs. Our theory reproduces the dependence of the excitation energies 
on the filling factors as measured in |14| reasonably well when v < 1/5. We also find that 
the shear modulus of the CF lattices goes down as the filling factor v = 1/5 is approached 
from below. This behavior is consistent with the experimentally observed increase of the 
threshold voltage for filling factors v — > 1/5 p3f (if one interprets the results in terms of 
"weak pinning" 

For v > 1/5 we will show in section V that the energy landscape becomes very flat, with 
many different lattice structures becoming nearly degenerate in energy. Not coincidentally, 
the convergence of our Hartree-Fock procedure is very poor in this region, and we are unable 
to identify the proper ground state in the clean limit. We have therefore presented two values 
for the gaps in this region, the upper one being for the triangular lattice, and the lower one 
being for a more oblique lattice. Neither gap follows the experimentally observed non- 
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monotonic dependence E g (u). (Note, however, the different slope below and above v — 1/5 
in Fig. 1). We believe the main reason for this is the following: Since there are many local 
minima with different lattice structures which are very close in energy, disorder may play 
an important role in real samples. The experimental gaps may also be dominated in this 
region by disorder effects. Apart from this one region of discrepancy, our numbers for the 
gaps are in reasonable agreement with experiments. 

The outline of the paper is as follows. In Section II we introduce the Hamiltonian 
formalism and show how the wave function (§) emerges naturally in this approach. In Section 
III we derive the expression of the electronic density operator in the CF representation and 
also calculate the Hall conductance for the CF lattice in the clean limit. In Section IV we 
formulate the HF theory using the electronic density derived in Section III. In section V we 
present the results and discuss their physical import. Some details of the calculations are 
relegated to the three appendices. 



II. SETTING UP THE HAMILTONIAN FORMALISM 



In this section we follow closely the exposition and notation of Murthy and Shankar |25 
Since most of the details are similar, we will mainly highlight the differences. We will follow 
Lopez and Fradkin |J in assuming that a good starting point for obtaining a perturbative 
solution of the N interacting electron problem moving in a uniform magnetic field B = —z B 
is the no n- interacting CS particle Hamiltonian (h — c — 1, m is the mass of the electron, e 
its charge) 

N , 

H cs = E ^H V * + eA *( r *) + a ^(r,)) 2 . (4) 
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The CS particle is obtained from the electron by attaching / flux quanta to it. This flux 
attachment is the origin of the Chern-Simons gauge field acs(rj), which is defined by 

Z^^£i> = £y (r _ rj) _ n(r) , (5 ) 
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where n(r) is the expectation value of the electron density in the ground state (CS particle 
and electron densities are the same). Unlike the incompressible FQH state (liquid state) 
for which the electrons have a uniform density, we are concerned with the case where the 
density depends on the position. The vector potential A*(rj) corresponds to the difference 
of the external magnetic field and the average field created by the attached flux tubes, or 
B*(r) = B - 2nln(r)/e. 

Our main assumption throughout this work is that for the appropriate field strength and 
electron density the electrons organize themselves into a density wave state such that the 
ground state expectation value of the density operator is a periodic function 

n(r) = n + £ 5n(G)e iG ' r , (6) 

G^O 

where G are the reciprocal lattice vectors. We also assume that the uniform component 
of the average density is much larger than any of the finite-G modulations, i.e. n ^> 5n. 
Effectively Sn/n is a small parameter in our theory. The effects of disorder are ignored in 
what follows. 

The complicated part of the CS Hamiltonian in (0) is the gauge field &cs( r i)- To get rid 
of it, one enlarges the Hilbert space Q by introducing a canonical pair of fields a(q), P(q) 
for every q 

[a(q),P(q')] = z(27r)V(q + q') (7) 

where q < Q = yAwn. Instead of working with the CS Hamiltonian, we introduce an 
equivalent Hamiltonian 

N 1 

H = £ ^H V * + eA *( r *) + a ( r *) + a cs(r*)) 2 , (8) 

where a(r^) = — iz x qa(rj) is a transverse vector field. We also define a longitudinal vector 
field P(r.j) = zqP(rj) (q is a unit vector in the q direction). This problem is equivalent 
to the original one provided we restrict our attention to states that are annihilated by the 
constraints x(q) = a(q) (q < Q). We will call states that are annihilated by these constraints 
physical states, that is, 
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X {<l)\% h ys>=0. (9) 

We will continue to use notation f° r the constraint operator in different representations. 

Using the fields a(q),P(q), a unitary transformation is then constructed which shifts a 
to absorb the Fourier components of SLcs{ r i) for q < Q. In the new representation, whose 
particles we call composite particles (CP), the Hamiltonian is 

TV , 

H cp = E + eA*(r<) + afo) + 27rZP( ri )) 2 . (10) 

We neglected all the q > Q Fourier components of the gauge field in deriving the above 
Hamiltonian, implying that our theory will not describe the motion correctly for large 
momenta (or short distances). In the following we will also be using the Random Phase 
Approximation (RPA) generalized for the case of the inhomogeneous densities 



j 

The constraints in the CP representation are given by 



£ e ik rj ^ n(2n) 2 5 2 (k) + 5n{-k). (11) 



x(q) = -^+p(q)-n(q), (12) 

where p(q) = J2f e _iq rj is the CP density operator and coincides with the corresponding 
electron operator. 

Before proceeding with any further transformations on the Hamiltonian (0), we will 
show how the trial wave function used to compute the energy of the correlated WC in 
emerges naturally within this approach |25||. The crudest approximation for the CP 



Hamiltonian is 

N i 

tf CP ^£ (_, Vl + eA*( ri )) 2 
2m 

Q 

71 _ _ 

+ ^ 5>(-q)a(q) + (27r/) 2 P(-q)P(q)). (13) 

In this expression we have assumed that the CP ground state can be regarded as largely 
homogeneous, so that the Sn is neglected in the definitions of the RPA (|II]) and vector po- 
tential A*(rj), which corresponds now to the uniform magnetic field. We have also neglected 



the coupling between the CPs and the oscillator fields a(q), -P(q). Since this Hamiltonian 
(|T3|) has been artificially made separable into a sum of the particle and the oscillator terms 
we can write down the ground state as a product wave function. The particles are moving 
in a uniform magnetic field B* = B — 2nln/e, so there is a degeneracy in this problem, but 
we assume that their ground state is crystalline 

^(M)=n<M r *), ( i4 ) 

i 

where 0^ (rj) are Gaussians centered on the lattice sites Rj similar to @, except that 
instead of magnetic length lo there is a new magnetic length Zq = (eB*)^ 1 ^ 2 . CF stands for 
composite fermions. The oscillator term describes N independent harmonic oscillators with 
the ground state 

Q 



^o SC ({q}) = n^ a(q)/4 ^- (15) 
q 



Using the constraints (|T^) we can eliminate the field degrees of freedom a(q) in favor of the 
Yi in the expression for the oscillator wave function QT5|) . Since the calculation is described 
in in great detail, we only give the final result for the projected oscillator wave function 

* osc ({ri})( a ( q )=27rZp(q)/ g ) = II l^i — e~ ^ Mzj |2/4i< > . (16) 

Here Zi — Xi + iyi is the complex coordinate and v = 2nn/eB is the filling factor. The 
approximate CS wave function is a product of (|T4| ) and (|1^) and agrees with the equation 



(2.12) given in Ref. [|0 



However, the fact remains that in this representation the oscillators and the particles 
remain strongly coupled. The oscillators are identified with the magnetoplasmons, which 
are high-energy degrees of freedom, while the particles will turn out to form the low-energy 
sector. Our next task will be to construct a canonical transformation to decouple the two 
sectors, so that we are left with a purely low-energy theory. 
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III. DECOUPLING AND THE ELECTRON DENSITY OPERATOR IN THE 

FINAL REPRESENTATION 



Before we turn to the technical details of the decoupling transformation, it is worthwhile 
to articulate the philosophy of our approach. If one were able to find the exact canoni- 
cal transformation, and implement it exactly, then one would be left with a final theory 
in which the fermions are purely low-energy objects, and the oscillators are purely high- 
energy objects. In particular, the oscillators should obey Kohn's theorem while all 
reference to the bare mass should have disappeared from the low-energy fermionic part of 
the Hamiltonian. In other words, the bare kinetic energy of the CFs should be quenched 
in the final representation. Finally, the projected electronic density when expressed in the 
final representation should obey the magnetic translation algebra appropriate to the LLL 
|29fl . Unfortunately, this program cannot be implemented fully in practice. What can be 
implemented is a sequence of transformations that achieves some measure of the above at 
long distance scales (small q). We will see that the oscillators do end up obeying Kohn's 
theorem, since this is a small-g property. Similarly, the magnetic translation algebra will be 
seen to occur in its small-g form. However, while the tendency for the quenching of the mass 
will be manifest, the mass depends on all distance scales, and its quenching cannot be shown 
within a long-distance approximation. Our approach will be to assume the exact quenching 
of bare mass, since we know this to be true in the LLL, and write the final Hamiltonian 
in the low-energy sector as a pure interaction term. Thus, while the proximate goal of the 
canonical transformation is to decouple the high- and low-energy parts, its ultimate goal is 
to obtain the electronic density operator in the final representation. 

We return now to the CP Hamiltonian flT0|) and the set of CP constraints (|12|) . Using the 
RPA approximation as given by (|TTD, the CP Hamiltonian can be recast into the following 
form 

N tt2 f—j Q 

^P = E^ + ^E(Aq) C f (q)+H 

i q 
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+-EE((*,» + MG)) 

q G 



X 



(At(q)A(q - G)$_(g - + h.c.) . (17) 



The first term in Eq. ([T7p is the CP kinetic energy, with FT,- = — zVj + eA*(i\,), which only 
depends on the particle degrees of freedom. The second is the coupling between the particle 
and the auxiliary field (oscillator) degrees of freedom, with c(q) = g_ J2j nj + e~ jq ' rj and the 
'destruction' operator A(q) = (a(q) + z27rZP(q)) / V 4nl p6| . The last term describes the 
oscillators and does not depend on the particle degrees of freedom. In order to decouple 
the high energy oscillators from the low energy CP's we need to compute the commutators 
between the newly introduced operators. It is straightforward to deduce from (|7|) that 

A(q),At(q')] =(27r) 2 5 2 (q-q'). (18) 

The commutator for the operator c(q) is found from the commutator of the canonical mo- 
menta Hz and Uy and then using the RPA approximation ([TT]) . The result is 

c(q), c^(q') ~ (— 47r/n5n(q — q') 
+ 2eB* (n(27r) 2 5 2 (q - q) + 5n(q_ - q'))) . (19) 

According to our assumption the average density is a periodic function, as in @, therefore 
the right-hand side of Eq. ( |19|) differs from zero only if the difference q — q' in ( |19"D is equal 
to a reciprocal lattice vector G. 

Our task is to decouple the Hamiltonian (|T7| ) by eliminating the term that couples the 
particle and the oscillator degrees of freedom. We will show how one can construct a canon- 
ical transformation that accomplishes this decoupling. Once the canonical transformation 
is found we can derive the electron density operator in the 'Final' representation (FR). Op- 
erators in the FR are expressed in terms of the CF coordinates. We will show below that 
there are good reasons to believe that Composite Fermions are the true quasiparticles in the 
FQH regime and the FR density operator represents the physical charge density. 
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The calculation is a straightforward extension of the procedure given in reference p5[ for 



the case of the homogeneous liquid, and we relegate the details of this calculation to Ap- 
pendix I. In what follows only the results of applying the transformation on the Hamiltonian 
(|TT|), the density operator and the set of constraints (|T2| ) are presented. 

By construction, the term coupling the particles and the oscillators is not present in the 
FR Hamiltonian. Substituting Equations (|66|) , ( |67|) and ( |70|) into the expression for the CF 
Hamiltonian ( |17|) we find that the oscillator term in the FR is equal to WcY^q A j! (q)A(q) 
with uj c = eB/m exactly as in the liquid, to order (Sn/n) 2 . This is a physically correct result 



because according to the Kohn's theorem |28| the limit co> c (q — > 0) should not depend on the 



electron interactions in the lowest Landau level. The kinetic energy in the FR is 
m A IL_IL+ A feB* irl , A 

t = y -i — 3 -± + y for,- 

y 2m y\2m m V j! j 
' Ec t (q)c(q) 



2mn 

+2^i E E ct (<l)c(q - G)MG)i( 9 -G) + . (20) 

PI G 

In an ideal calculation the particle kinetic energy should disappear in the FR; the electronic 
kinetic energy is subsumed into the oscillator term. As has been stated above, it is impossible 
to show this in a small-g approximation such as the one we are using. 

The electron density operator in the FR is obtained by solving the flow equation that is 
derived in a way that follows closely the calculation for the kinetic energy T leading to the 
Eq. fl69|) . The result of the integration of the flow equation is 

p(q, A) = p(q) + p (q, A) + Q £ U(q - G)g_(g^G) + + h.c. 

W7r/(l + /i 2 ) g N V 



where /i 2 = 1//^ — 1. The FR operator Po(qA) is the leading term in the perturbation 
expansion of the density in the parameter 8n/n and is formally identical p7| to the density 



operator in the case when the average electron density is uniform |25 
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2^/(1 + /i 2 ) 

V c + h c (22) 

47rZr^y(l + /i 2 ) 

It is now straightforward to get the FR expression for the set of constraints (|12|), just 
by using the previously determined FR operators A(q, A) (see Appendix I, Eq. fl66|) ) and 
p(q, A) (Eqs (pl|) and (|22])). The expression that results is 

X (q, A) = Xo(q , A) - Mq) + 8 ;^; V 1 + ^ 

EMG) (c(q- G)g_(g^G) + + h.c.) . (23) 



x 

G 



Here again Xo(q> ^) is the P ar t of the constraint that corresponds to the case of the uniform 
average density and is given by 



Xo (q, A) = p(q) + 9(V ^^ 2 1} (c(q) + h.c.) . (24) 



The main observation about Eq. ( p3] ) is that oscillator degrees of freedom cancel out up 
to order 5n/n, implying that the constraint acts only on particles. This reassures us of the 
self-consistency of the decoupling scheme, since there is no use decoupling the high and low 
energy modes in the Hamiltonian if the constraint still non-trivially couples them. 

Because the particles are confined entirely to the lowest Landau level, one expects the 
physical charge density operator to obey the magnetic translation algebra. However this 
is not true for the density operator defined by Eq. (21). We are allowed to modify the 



definition of p(q, A) in CP representation by adding to it any multiple of the constraint, 
since in an exact calculation in the physical states the constraint is equal to zero. Following 
the same approach as in the liquid states |25| we try the linear combination 



p(q)-l/(^ + lMq). (25) 

This operator has the virtue that its FR matrix elements are of order q 2 or higher, consistent 
with Kohn's theorem []2S|] . The FR expression of the density operator fl2"5] ) (which we will 
call the preferred density) is 
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P(q) = ^TP(q) " ; , ,? . * (c(q) + h-c.) + Mq 



/i 2 + l rv ^ 4vr/n(l + /i 2 ) v ^ ' v /i 2 + l 
g (l-vT+^) 



£ Sn(G) (c(q - G)g_(g - G) + + h.c.) . (26) 



47rln(l + /x 2 )/i 2 

The calculation of the commutator of the preferred density operators to first order in 
8n/n gives 

[p(q), p(q')] = i*o(q x q')p(q + q') 

+«2(q x E ^^Xo(q + q' - G, A). (27) 

Here Xo(q ; A) is the constraint to zeroth order (in 5n/n), Eq. (f24j). We conclude that the 
magnetic algebra is satisfied for physical states that are destroyed by the constraint. This 
is a weaker result than was obtained in the translationally invariant case, but nevertheless 
still preserves the equivalence of this theory to the original electronic theory in the LLL at 
long distances. 

The preferred density encodes many nonperturbative features that are known to be true 
for the original electronic problem. It shows the correct fractional charge of the quasiparticles 
0, obeys the magnetic translation algebra in the small q limit P5| , and has matrix elements 



of order q 2 or higher from the ground state. Thus it is plausible that simple HF calculations 
with this density will capture the essential physics in the FQH regime. Certainly, this 
expectation has been borne out in calculations for the liquid states ||. 

Thus all the features of the translationally invariant Hamiltonian theory, namely compli- 
ance with Kohn's theorem, simultaneous decoupling of the Hamiltonian and the constraints, 
and a preferred density that obeys the algebra of magnetic translations, carry over for the 
non-uniform electronic density state. 

It is of interest to determine the Hall conductance a xy for the non-uniform average density 
state. In the clean limit, when there is no external potential, translation invariance implies 
that a xy = ue 2 /h. Our theory does indeed predict this in the clean limit. In order to see 
this we need the FR expression for the current. Starting with the electron current, and 
eliminating the CS vector potential we find the following CP current 
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1 2m/7r/ 2\/7r/ 

(Jcp) + (q) = - E n i+ e-«™ + -^A(q)g + + £ 5n(G)A(q - G)(g - (28) 

The first term in Eq. ([28|) is just a definition of the operator c(q) so to get the FR expression 
for the current we have to substitute the FR operators A(q, A) and c(q, A). The FR current 
consists of two terms, one that does not contain factors of 5n and another proportional to 
Sn/n. The first term is identical to the FR current for the uniform average density state 
and as shown in [25|| depends only on the oscillator degrees of freedom. To calculate the 
first order contribution to the current we use expressions (|B6D and (^?p that give the first 
order corrections in 6n/n for the operators A(q, A) and c(q, A) respectively. We find that it 
is also independent of the CF coordinates. Both terms add up to 



2nJirl(l + a 2 ) Jvr/M+u 2 ! — 

■4(q) = — ^ —Am + + ^ — E *n(G)A(q - G)(g - G) + . (29) 

m m g 

Because the current in the FR depends only on the operators that represent the oscillators, 



we can ignore the particle sector in the conductance calculation. The argument in Ref. [25 



for the uniform average density case goes through and gives a xy = ve 2 jh in the limit q — > 0, 
which is the correct unquantized Hall conductance in the clean limit. Note, however, in the 
presence of disorder it is believed that a xy — > 0, a xx — > such that p xy is its classical value. 
A complete theory including disorder effects is currently nonexistent, and we will confine 
ourselves to the clean limit in the sequel. 



IV. HARTREE-FOCK APPROXIMATION 

Having determined the correct canonical transformation by decoupling the CS Hamilto- 
nian, we shift our focus to the Coulomb interaction that was hitherto ignored. At the small 
filling factors for which the Wigner Crystal occurs, the lowest Landau level approximation 
is appropriate. Now the LLL electronic Hamiltonian is given by 

h = ^ E (V (q)p(q)p(-q) - v( q )e~^ 2 P (o)) . (30) 
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Here S is the area of the system and the second term is a result of un-normal-ordering the 
original electronic Hamiltonian. This is necessary since we need to have the full operator p 
in order to write it in the final representation. The two-dimensional Fourier transform of 
the Coulomb potential is suppressed at large momenta q by multiplying it with a Gaussian 

V(q) = — e"^, (31) 

where a parameter A may be used to interpolate between the pure Coulomb potential and 
the Coulomb potential that is effective only when the the distance between two particles 
is larger than A. e is the dielectric constant. We emphasize that the form of the Fourier 
transform of the potential given in Eq. ([H]) does not accurately describe the effect of the 
sample thickness, but is rather chosen for illustrative purposes, since it is computationally 
convenient. The second term in Eq. (^) is a constant that will be ignored in what follows. 
The first term of the electronic Coulomb interaction when transformed to the CF coordinates 
will serve as our model Hamiltonian 

H = ^J2V(q)p(q)p(-ci). (32) 

The density operator p(q) is given by Eq. (p6|) . It is useful at this point to rewrite it so 



that the dependence on the modulated average density 8n is explicit |30 



p(q) = (1 - c 2 )p(q) - il 2 E q >< "/ '"• r - 



2 MG)q- G i(G _ q) . r . 

N G 2 



-cr E 



+ iMl^MG)^ qxn , e , (G - q ,. rj (33) 

Here II* is the momentum operator that corresponds to the uniform average density case 
(IT* = Uj (8n = 0)) and c = y/lu. It will also be convenient to have separate symbols for 
the different orders of the 5n in Eq. (|33|), so we write p(q) = po(q) + J2g ^ n (G)pi(q, G). 

The Hamiltonian (^) describes a many-body CF problem that we will treat within the 
Hartree-Fock approximation. We justify the use of the HF approximation by arguing that 
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the CF is the true quasiparticle in the FQH regime. It will be assumed throughout this 
study that the average density modulation is small compared to the uniform background 
(Sn/n <C 1) so as a convenient basis we will choose the wave-functions of the free CF moving 
in the uniform magnetic field B* = B — 2-Kln/e. The Landau gauge will be used in what 
follows. Wavefunctions will be denoted as \n, X), with n as a CF Landau level index and X 
as a kinetic momentum component in the y direction. 

We will now derive the HF Hamiltonian. The model Hamiltonian (|32l ) in the \n, X) basis 
may be written 

ZD q mXi,...,n 3 x 3 



t 

71] 

5>(q) £ (nAl^q) (714X4) 



x (n 2 X 2 | p(-q) |n 3 X 3 ) e^c^,^ 
1 

+ 2S 



Q niJfi ,...,124X4 

x (n 2 X 2 | p(-q) |n 3 X 3 ) 

C ni ,Xi C ri2 ,X 2 C «-3^3 C n 4 X 4 j 

(34) 

where c^ nX (c n; x) is the CF creation (destruction) operator. The usual HF pairings are made 
in the two-body term of the Hamiltonian (|3"4]), giving two contributions - a direct and an 
exchange term. Because the Xj dependence of the density matrix elements in Eq. ( [34]) is 
very simple 

P (q) |n 2 X 2 > = (m| p (q) \n 2 ) e-^ +x ^ 2 5 Xl ^ 2 _ gylf (35) 

(mXi| pi(q, G) |n 2 X 2 ) = (m| Pi(q, G) |n 2 ) 

Xe -i(,- Gl )(x 1+ , 5 )/ 2 ^ a2 _ fe _ Gj)i;j (36) 

(here Zq = (eB*)~ l l 2 is the CF magnetic length) one can eliminate all the dependence on 
the momentum Xi in the Eq. (0). To this end we introduce the operator 

A»'(q) = ^£e-^^^^'5 2 / 2 ct iX c n , x+g!/i52 , (37) 
9 x 

where g is the degeneracy of a Landau level. A nn /(G) is the order parameter of the density 
modulation corresponding to the wave- vector G. Note that A nn /(q) has the property 
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£A nn (0) =u. 

n 

After doing sums over Xi we find the following contributions to the HF Hamiltonian 

1. A one-body term, zeroth order in Sn: 

<4^ 7(q) S K|po(q)|n 2 )(n 2 |p (-q)K)A„ in3 (0). 

q n\,ni,nz 

2. A one-body term, first order in 5n: 

C = ^E^q) E «ni|^(-q,G)|n a >(n 2 |pb(q)|n 3 > 

D q,G ni,n 2 ,n 3 

+ (nil Po(-q) |n 2 ) (n 2 | pi(q, G) |n 3 » 
x5n(G)e^ xG / 2 A nin3 (-G). 

3. A one-body term, second order in 5n: 

H ll = {o E V(q) E (ni|pi(-q,G)|n 2 >(n 2 |p 1 (q,G 1 )|n3> 

q,G,Gi n\,ni,nz 

x<5n(G)5n(G 1 )e^ 2 ( qx ( G+Gl )- GxGl )/ 2 A nin3 (-G - d). 

4. A two-body term, zeroth order in 5n, direct and exchange contributions: 

^° = ^E^(G) E KIPo(-G)K)(n 2 |p (G)|n 3 ) 

G ni,...,ri4 

x(A nin4 (-G))A„ 2 „ 3 (G)-E^(q) E KIPo(-q)K) 

q,G ni,...,n 4 

x| (n 2 |p (q) |ns> (A ftl „ 3 (G)> A n2n4 (-G)e^ Gx <>. 

5. A two-body term, first order in 5n, direct contributions: 

<, = ^E E (^(-G - d) (ni| pi(G + Gi, G) \n 4 ) 

G,Gi ni,...,U4 

x (n 2 \ p (-G - d) |n 3 ) + V(-G0 (m| jo (G 1 ) |n 4 ) 

x (n 2 | pi(-Gi, G) |n 3 )) tfn(G) (A nin4 (G0) A nan3 (-G - Gi). 
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6. A two-body term, first order in Sn, exchange contributions: 

< e = -f E E nq)(Klpi(-q-G)|n 4 >{n 2 |po(q)h> 

q,G,Gi »l »M 

Xe « 5 =»(Gx(q-G 1 )/2+G lX q) + ( ni |p (_ q ) | W4 ) 

X (n2|i9i(q,G) | n3 )e^ 2 (Gx(q + G 1 )/2 + G 1 x q) ^ 

x5n(G) (A nm3 (G 1 )) A„ 2 „ 4 (-G - d). (44) 

7. A two-body term, second order in Sn, direct contribution: 

h]L = 9 ^ E E n-G-co 

° G,Gi,G 2 ni,...,n 4 

x (m| p x (G + Gi, Gi) |n 4 ) (ra 2 | pi(-G - G x , G 2 ) |n 3 ) 
x5n(G 1 )5n(G 2 ) (A ni „ 4 (G)) A„ 2 „ 3 (-G - G x - G 2 ). (45) 



A two-body term, second order in Sn, exchange contribution: 



HlL = ~ E V(q) E (nxiM-q.GOM 

q,G,Gl,G2 m,...,Ti4 



x (n 2 | pi(q, G 2 ) |n 3 ) (A nin3 (G)) 5n(G 1 )5n(G 2 )A n2 „ 4 (-G - Gi - G 2 ) 

xe ii* 2 ((Gi+G 2 )xq+Gx(Gi-G 2 )-G 1 xG 2 +2Gxq)/2 



The matrix elements of the density operator can be calculated using the formulas (|7T|), ([72[) 
and (^) that are given in the Appendix II. The momenta Gj run over a discrete set of 
reciprocal lattice vectors. The momentum q is a continuous variable. The summation over 
q in those terms of the Hamiltonian where it appears can be done in a closed form as we 
show in Appendix III, because the potential in Eq. (|3l"l) was chosen so that these integrals 
could be performed analytically. 

We will group all the non-operator entries in the Equations (f39j)-(f4"6l) under the notation 
U nin , 2 (G) (renaming the dummy summation variables where necessary) and represent the 
HF Hamiltonian in a form convenient for further discussion, 
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Hhf = 9 E *WG)A nina (G). (47) 

G,ni,n2 

Obviously U niJl2 (G) depends on the expectation value of the order parameter operator 
A mn2 (G) both directly and through the density modulation Sn, because 

5n(G) = g E (m\ p (G) \n 2 ) (A mn2 (G)> 

m,n2 

+2 E (n 1 |p 1 (G,G 1 )|n 2 )5n(G 1 )(A nin2 (G-G 1 )}. (48) 

Gi,ni,na 

In the z'th iteration of the numerical HF procedure A Tlin2 (G) is calculated using the solution 
of the z — l'st (previous) iteration. The density modulation is then calculated as a numerical 
solution the system of the linear equations defined in Eq. ([4"8"[). 

Having found the HF Hamiltonian, we can solve it to find the single-particle spectrum of 
the many-body system. We will assume that the CF form a Wigner lattice with one particle 
per unit cell. The reciprocal lattice constant for a triangular lattice is given by 



However, we will find that in some regions of filling factor the triangular lattice is not the 
ground state, and we will explore other lattice structures. 

We have used two different schemes to perform the calculation, one due to Cote and 



MacDonald [31| and the other due to Yoshioka and Lee [32]. Below we will outline the 
essence of each of these methods. 

The method by Cote and MacDonald (CM) starts from the single-particle Green's func- 
tion which they define as 

G nin2 (X u X 2 ,r) = - (T Clll , Xl (r)i iX2 (0)) , (50) 

here T is the time-ordering operator. The relationship of the Green's function Fourier 
transform to the physically relevant expectation value of the order parameter is 

(A nina (G))=G nam (G J r = 0-) 
= \ E G n2ni (X 2 ,X l ,r = 0-)e-* G ^ +x ^ 2 5 XuX2 _ Gylf . (51) 
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The set of crystal order parameters (A nin2 (G)) is then used to find the ground-state energy 

Ehf = - E tWG)(A„ in2 (G)>, (52) 

G,ni,ri2 



where e\ = 1 if it multiplies those terms of U nin2 (G) that are given by the Equations (|39|)- 
(one-body terms) and ei = 1/2 if it multiplies terms that are given by the Equations 
(f£2|)-(|46|) (two-body terms). The excitation gap E g , also called activation energy, can be 
deduced from the chemical energy and the single-energy density of states d(E) which is 
related to the Green's function through 

d(E) = — V %G nn (G = 0, iuj -^E + iS), (53) 
^ n 

here 5sG nn is the imaginary part of the operator G nn , Uj are the Matsubara frequencies, 5 
is a small smoothing parameter. 

All of the above is predicated on knowing the Green's function. We derive the Green's 
function equation of motion in the usual way by taking the commutator of the Hamiltonian 



( fj7|) with a single particle destruction operator c n x 



iojj + j^j G nin2 (G,Uj) 
~ E t^WGi - G)G n3n2 (Gi, ^)e iGxGl ' 52 = 5 niin2 8 G> o, (54) 

G,n 3 n 

where \x is the the chemical potential. The system of Equations ( |54l) is solved for the Green's 
function by diagonalizing its left-hand side with respect to the indices n 3 and Gi. One can 
find the expectation value of the order parameter and the density of states once the chemical 
potential is known. The chemical potential in turn is calculated by filling up the correct 
number of states, that is, by using Eq. (p8|). 

The numerical iterative scheme starts by assuming a Gaussian form for the order pa- 
rameters. (The exact expression depends on the filling factor and the state that is being 
constructed and will be discussed later). This initial set of (A nin2 (G)) is then used to com- 
pute the effective potential U niri2 (G). Next the equation of motion is solved to get a new set 
of order parameters and the process is repeated until the (A nx „ 2 (G)) converge with some 
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prescribed accuracy. Another way to check the accuracy of the numerical solution is by 



using the following useful sum rule BT] 



E (A nin2 (G)) 2 = v. (55) 

G,ni,ri2 

The idea behind the second method, that of Yoshioka and Lee (YL), is to diagonal- 
ize the one-body HF Hamiltonian that can be rewritten in terms of the CF creation and 
destruction operators as 

H = E U nin 2 ( G ) e ~ tGxXc i u X-G y l* 2 /2 C n2,X+G y l*y2- (56) 
G,ni,n 2 ,X 

We assume that the CF form a periodic lattice with primitive translation vectors of the 
reciprocal lattice that are given by Qi = (Qo,0) and Q2 = Qo(p/q,a). The first unitary 
transformation on the Hamiltonian fl56|) is defined by 

* ni ,x,Y = ^=E e- lsaQoY c nuX+saQolf , (57) 

V b m s= q 

here s m = L/aQol^, L is the linear dimension of the system, < X < aQol^ and < Y < 
2ii/aQo. After making the transformation fl5"T|) the Hamiltonian is 

tt — ST TT (C^\„~iG x X+iG y Y+iG x G y ll 2 /2 n ] /roN 
a - 2^ U niri2 {G)e vo ' a ni , X ,Y a n 2 ,X,Y+G x l* 2 - (5°) 

G,m,n 2 ,X,Y 

The variable Y in Eq. (|58| ) is coupled through G x Iq 2 . If this number is commensurate with 
27r /aQo, which is the period of the variable Y, then we can simplify the Hamiltonian even 
further. Suppose then that the parameters are such that NQol^/q = M2tt /aQo, with M 
and N integers. We then introduce a new operator 

K,j,X,Y = a n,X,Y+jQ lf/qi (59) 

where 1 < j < N and < Y < l^Qo/qM. After inserting Eq. (^) into ([58]) the expression 
for the Hamiltonian is 

# = EE E U nin2 (G)e-^ x+ ^ Y 

G X,Y m,j,n 2 ,k 

v p iQoQyjlo 2 /q+iG x G y l* 2 /2i\ i r / fir) \ 

Xe °n 1 ,j ) X,Y °n 2 ,k,X,Y°k,j+Q x q/Q - {VU) 
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For every pair (X, Y) that takes values in the rectangular domain defined earlier the Hamil- 
tonian flBTf ) can be diagonalized in the indices and ri2,k. The single-particle energies 
E ni j(X, Y) that result thereby are continuous in the variables X, Y and form energy bands. 
There are n m N energy bands (where 1 < ni, n% < n m ) and there is a large energy gap be- 
tween the lower Mg'th and Mq + l'st bands. The CF state that has the chemical potential 
in this large gap should have the lowest energy. Such a state occurs when the CF filling 
factor is v* = qM/N, where we have defined the CF filling factor as u* = I^V^n. It is easy 
to derive the expression for the expectation value of the order parameter, by applying the 



transformations (p7|) and (|59|) to the definition to get 



L X,Y j,k 

xe *QoQ y] if/ q+ >G x G y if/2 ^ ^ 5kd+Qxq/Qo . (61) 

As in the case of the previous method we find a solution to the HF problem by iterating 



until the order parameters converge. We calculate the ground state energy using Eq. ( p2|) 



and the excitation gap E„ as a smallest separation between the Mg'th and Mq + l'st bands 



While the method of CM is numerically efficient it is sometimes difficult to extract the the 
excitation gap from the smoothed density of states. There is no uncertainty in determining 



E g when the method of YL is used. 



V. RESULTS 



Our experimental motivation is the work by Jiang et al. |T3|,[T4|| where the transport 



properties were measured around v = 1/5 Landau level filling. In |T3 an insulating phase 
was identified just above v = 1/5 at v — 0.21 by observing a large peak of the longitudinal 
resistance R xx as a function of the external magnetic field. The activation energy was 
estimated from the Arrhenius plot at E g ~ 0.63 K (with B w 20 T). The striking thing is 
that the magnitude of the activation gap compares very poorly with the results obtained 
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from HF for the usual electron solid. The excitation energies for the triangular electron 
lattice with one particle per unit cell are given in Table 1. We use the modified Coulomb 
potential given by ( |3~TD and present results for different values of the thickness parameter 
A. The calculation was done in the lowest Landau level approximation and for A = it 
reproduces previous results J3"2"| . The energies are given in units of e 2 /el . In the same units 
the experimental result is E g ~ 2.8 x 10~ 3 e 2 /e/o, at least two orders of magnitude smaller 
than the theory. 

Table 1 Electron WC ground state and activation energies for different values of A. 



A 





Z /2 


k 


3l Q /2 


Ehf 


-0.3220 


-0.3137 


-0.2859 


-0.2413 


E g (e 2 /d ) 


0.4728 


0.5080 


0.5080 


0.4468 



We expect some reduction in the value of E„ when the relaxation of the lattice is ac- 



counted for |35| , |34]| , but it is difficult to believe that this correction would nearly exactly 
cancel the unrelaxed excitation energy. Besides one would not expect the E g {y) for the 



electron WC to be non-monotonic as observed experimentally [14 



Now we proceed to carry out our program of considering crystals of CFs. 



A. Crystals of Composite Fermions with two vortices attached 



Let us examine how well CFs with / = 2 describe the experimental situation. An elec- 
tronic filling factor of v corresponds to a CF filling factor v* = vj{\—lv) = 1/3. So the lowest 
CF Landau level is partially filled and it is reasonable to expect that the composite fermions 
form a lattice. As in the electron solid calculation only the lowest CF Landau level is kept 
(n m = 1). Keeping two CF Landau levels (n m = 2) we find similar results, indicating that 
including more Landau levels does not influence the calculation. Because the CF and elec- 



tron effective potential Uoo(<l) in Eq. (47) are different momentum functions it is not obvious 
that the CF lattice is triangular as is the case for the electron lattice. The functional form 
of Uoo(q) may be suggestive in that respect. One expects it to have a minimum at about the 
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momentum q equal to the shortest reciprocal vector. This argument cannot be exact in our 
theory because Uqq(G) depends on 5n and is reevaluated self-consistently in every iteration. 
However, since the density modulations are small, we expect that a good approximation 
to f^oo(G) can be obtained by keeping only the 5n- independent terms given by Equations 
and ( f4"2"D in the HF Hamiltonian (the term given by Eq. ( |3"9"D is a constant). In that 
case the approximate effective potential can be expressed as [/00(G) = Wo(G) (A o(G)), 
defining the effective interaction Wq{G). We display the plot of this effective interaction in 
Fig. 2 for different values of the parameter A. Whereas for the v* = 1/3 triangular lattice 
we expect a minimum at about |q| /g ~ 1-56, the minimum for CF effective potential is at 
much smaller wave-vector, more so for a small A. This is why we do not limit ourselves to 
the triangular lattice but calculate the ground state energies along with the E^'s for three 
oblique (including triangular) lattices. The results together with the primitive reciprocal 
lattice vectors bi, b2 are given in Table 2. Every lattice is rescaled by an overall factor that 
makes the volume of the unit cell a constant equal to lul^jv. 

Table 2 / = 2 CF lattice ground state and activation energies for different values of A 
and different unit cells. 





A 





l /2 


lo 


3/o/2 


bi = (1,0) 
b 2 = (0.5, y/3/2) 


Ehf 


-0.33 


-0.34 


-0.30 


-0.24 


E g (e 2 /d ) 


0.08 


0.13 


0.13 


0.14 


bi = (l,0) 
b 2 = (0.5, x/3) 


Ehf 


-0.36 


-0.33 


-0.29 


-0.24 


E g (e 2 /d ) 


0.06 


0.07 


0.12 


0.12 


bi = (1,0) 
b 2 = (0.5,3^3/2) 


Ehf 


-0.37 


-0.33 


-0.29 


-0.24 


E 9 (e 2 /d ) 


0.04 


0.05 


0.09 


0.10 



We find that for A = composite fermions prefer the elongated lattices to the triangular 
one. For larger values of A the triangular lattice is favored. The results for the excitation 
energy are somewhat closer to the experimental value but still too large. 

The disagreement between the theory and the experiment is not only in the magnitude 
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of the activation energy but also in its dependence on the filling factor. In our theory 
with I = 2 the function E g (i/*) varies slowly and is monotonic around the CF filling factor 
u* — 1/3 (y — 1/5). Fig. 3 gives this dependence for the triangular lattice with A = 31q/2. 
The experimental function (see Fig. 3 in JTJJ) has a sharp peak between the filling factors 
v ~ 0.22 and v ~ 0.21 and for v < 1/5 it rises sharply and saturates at lower filling factors. 
Let us turn to CFs with four flux quanta to see how the results compare with experiments. 



B. Crystals of Composite Fermions with four vortices attached 

The behavior of the experimental gap with v fits in more naturally within the CF model 
with I = 4. When the filling factor v < 1/5 the lowest CF Landau level is being populated 



and a CF quasiparticle lattice is assumed to be the stable state fl36fl . On the other hand 
when v > 1/5 the second CF Landau level is being populated, and one naturally expects 
some difference in the behavior of the gap in the theory. We will see that this expectation 
is realized, but not in complete agreement with experiments. 

Numerical constraints limited our HF basis to the two lowest CF Landau levels (n m = 2). 
The initial seed used in the HF procedure that converged to the correlated WC was 

h ,* e -G 2 l* 2 /i if ni>n2 — 

otherwise. 

When v > 1/5, the second CF Landau level is partially filled. Again we assume that the 
CF lattice is formed so the initial seed that we use in this case is 



(A ni n 2 (G)) 



(A„ in2 (G)) 



e 



-G»tf»/4 if ni ,n 2 = 



y - l) e - G2/ (» 2 / 4 ifm,n 2 = 1 



otherwise. 

Our results for the activation energy are presented in Fig. 1. The value of the parameter A 
is 3Zo/2 (the results for A = 1$ are very similar) and we assume that the lattice is triangular. 
A magnetic field of B = 20 T was used to convert the energy units to Kelvin, in order to 



compare to the work of Jiang et al [ |l4j1 . 
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First we will discuss the results for v < 1/5. The excitation gaps that we obtain are 
generally comparable to the experimental values. We also reproduce a correct E g (y) depen- 
dence here (see the left half of the Fig. 1). However we do not observe saturation towards 
the lower filling factors. This may be an indication that perhaps CFs with I = 4 are not the 
quasiparticles at very low fillings. 

Another experimental probe supporting the crystalline nature of the insulating state 
is a / — V measurement ||23|| . Nonlinear I — V curves have a threshold voltage at which 



the differential resistance drops off that can be interpreted as a depinning of a weakly 
pinned Wigner Crystal [p4fl . As the filling factor is varied the threshold voltages increase 
approaching the FQH state at v = 1/5 both from above and below. This finding could 
be a consequence of a lattice getting less rigid as the FQH state is closer p4] , |37[| . We 
have calculated the shear modulus of the CF lattice for several fractions v < 1/5. We first 
compute the ground state energies of a triangular lattice with the primitive reciprocal lattice 
vectors bi = (1, 0), b 2 = (0.5, V3/2) and a deformed lattice such that it primitive reciprocal 
lattice vectors are bi = Q (1,0), b 2 = Qo(0.5, 3a/2/4) (oblique lattice), with Q chosen so 
that the area of the Brillouin zone is equal to that of the triangular lattice. Then the shear 
modulus fi is proportional to the difference of the ground state energies. The results are 
presented in Figure 4 for A = 1.5/q- We observe that the lattice is indeed becoming softer as 



v — > 1/5. This conclusion is consistent with the experimental results |23| interpreted using 
the collective pinning theory |pi,|37|. 



For v > 1/5 the gaps for the triangular lattice, while being in the same range as their 
experimental counterparts, do not show the correct dependence on v close to v = 1/5 (see 
the right half of the Fig. 1). We find that in this case the triangular lattice is not the lowest 
energy solution to the HF equations. Fig. 5 gives the HF energies of several lattices for 
a fraction v* = 6/5 (that corresponds to v — 0.206 . . .). The lattices that we consider are 
deformations of the triangular lattice obtained from it by changing the angle 9 between the 
reciprocal lattice vectors hi and b 2 such that |bi| = |b 2 | and the volume of the unit cell 
remains a constant. Fig. 6 gives the corresponding activation energies. As is apparent from 
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Fig. 5, the lattices with small angles 9 are more stable within our HF scheme. The smallest 
angle for which the iterations reliably converge is 9 = n/Q. The dependence of E g on v is 
presented for the triangular and the 9 = 7r/6 lattices in Fig. 1. 

At this point we have run into an intrinsic limitation of the Hamiltonian theory: Since 
the exact transformation between the electronic coordinates and the CF coordinates in not 
known, the Hamiltonian itself is not known exactly. This means that we should not take 
the ground state energies that are predicted by our theory too seriously. Note also that the 
differences in ground state energy between the different lattice structures are very tiny so any 
conclusion concerning the stability of one lattice compared to another should be taken with 
a grain of salt. We still can estimate the "shear modulus" for this class of lattices by looking 
at the difference in ground state energy between the triangular and square lattices. This 
leads to an estimate of /i ~ 2 x 10 -5 ^, an order of magnitude smaller than for v < 1/5. This 
means that the CF lattices are very soft for v just above 1/5, and disorder may potentially be 
very important in this case. As the filling factor increases the situation remains qualitatively 
similar but the differences in energy decrease. The HF energies for the triangular and square 
lattices are presented for several filling factors in Table 3. The activation gaps for these two 
lattices are essentially the same. 

Table 3 Comparison of I — 4 CF square and hexagonal lattice ground states for filling 
factors v > 1/5. Energy in units of -y-. 



V 


0.2069 


0.2105 


0.2143 


0.2174 


0.2195 


Ehf square 


-0.265046 


-0.265242 


-0.265716 


-0.266248 


-0.266659 


Ehf hexag. 


-0.265040 


-0.265232 


-0.265708 


-0.266247 


-0.266658 



VI. CONCLUSIONS, CAVEATS, AND OPEN QUESTIONS 

Two-dimensional electron gases in high magnetic fields offer the best conditions for the 
realization of the Wigner Crystal, since the magnetic field tends to localize the electrons. 
However, electronic correlations play a dominant role in the LLL because the kinetic energy 
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is degenerate. Attempts at describing the Wigner Crystal using uncorrected, or weakly 
correlated states of electrons ||12j| , do capture some of the essential physics, such as the filling 
factor at which the Laughlin liquid becomes unstable to the Wigner Crystal. However, these 
theories fail to capture the correct structure of the excitation spectrum, and predict gaps 
that are two orders of magnitude above experimental observations. 

Since Laughlin- Jastrow correlations are the essence of the fractional quantum Hall liquid 
states PJH, it is natural to hypothesize that they are important in the Wigner Crystal state 
as well. The first step in this direction was taken by Yi and Fertig who studied the 



ground state energy as more and more vortices were attached to electrons forming a Wigner 
Crystal. They found that indeed the Wigner Crystals with vortices had better energies than 
the uncorrelated or weakly correlated crystals ||20|| . 

Unfortunately, ground state energies cannot be probed directly in experiments. It is 
desirable to have predictions for observable physical properties that can distinguish between 
competing ground states. Calculating physical properties in a strongly correlated state is 
notoriously difficult. The Composite Fermion picture || achieves the miracle of transforming 
a strongly correlated electronic problem into a weakly correlated problem of CFs. In the 
years since the discovery of the FQHE, much progress has been made in developing field- 
theoretic schemes which have predictive power fH]]. The latest in this long line of approaches 
is the Hamiltonian formalism 0, which has had reasonable success in computing gaps, 
magnetoexciton dispersions, and finite temperature properties for the liquid states ||. 

In this paper we have partially accomplished the goal of computing the physical properties 
of a strongly correlated Wigner Crystal. Based on an extension of the Hamiltonian theory 
to account for the nonuniform density in the crystal state, we were able to compute gaps in 
the correlated crystal. 

Our results show that qualitatively and semi-quantitatively, a Wigner Crystal state of 
CFs with four flux quanta attached offers the best description of the phenomenology of the 
high-field Wigner Crystal near v — 1/5. In particular, our predictions for gaps are within 
a factor of 2 of the experiments in the entire regime of interest. Our predictions show a 
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different behavior for v < 1/5 and v > 1/5. While the theory has some discrepancies with 



the data for v just above 1/5, we believe we understand why this might be the case: 
Different lattice structures have very similar energies in this regime, and are very deformable. 
Consequently disorder is expected to play a dominant role in determining the configuration, 
and hence the gaps, in this region of v. Finally, we are able to estimate the shear modulus 
of the crystal above and below 1/5, and we find them to become softer as 1/5 is approached. 



This is consistent with the increase of the threshold voltage for nonlinear transport p3| , a 



standard feature of the theory of collective pinning p4| , p7 . 

Before we close, some caveats must be noted. An intrinsic limitation of the Hamiltonian 
theory is that the Hamiltonian is known only approximately. Thus the ground state 
energies are not to be taken too seriously. This implies that this theory does not offer a 
trustworthy way to find the lowest energy state. The strength of the Hamiltonian theory 
lies in the fact that if the nature of the state is known, the theory allows the computation 
of gaps, magnetoexcitons, and even finite temperature properties ||. With this in mind, let 
us note that we have not carried out an exhaustive search in the space of all possible states. 
We have confined ourselves to crystals with one CF per unit cell. While we did explore 
crystals other than triangular and square for v > 1/5, we kept the two primitive reciprocal 
lattice vectors equal in magnitude. It is possible that some other crystal state that we have 
not explored is the actual ground state in the clean limit. However, for the experimental 
observations this point is likely to be moot, because disorder probably plays a dominant role 
in this region of v. 

Many open questions remain. The most important, and the most intractable, is the 
influence of disorder. Disorder will cause lattice deformations, dislocations and other defects. 
In a crystal of CFs, density variations are expected to produce a corresponding variation in 
the effective magnetic field. Thus, a random potential leads indirectly to a random magnetic 
field. In principle, the formalism we have developed here to deal with nonuniform density 
could be generalized to incorporate disorder, but the implementation appears difficult. In 
particular, it is difficult to visualize how the nonperturbative effects of disorder (localization 
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of almost all states, changing <j xy from ve 2 jh to 0, etc.) would emerge in a straightforward 
manner. 

Another open question is the evolution of the Wigner Crystal state with temperature, 
which could be explored in the clean limit along the line of reasoning laid out in this work. 
In particular, it would be of interest to obtain a prediction for the transition temperature 
between the Wigner crystal and the (presumably liquid) high-temperature state. 

We hope to pursue these and other topics in future work. 



VII. ACKNOWLEDGEMENTS 

We are grateful to the NSF for partial support of this work under grants DMR-9870681 
(R.N. and H.A.F) and DMR-0071611 (R.N. and CM.), and the Center for Computational 
Sciences (R.N.) at the University of Kentucky. 



VIII. APPENDIX I 

In this appendix we will construct the canonical transformation that diagonalizes the 
Hamiltonian (|T7| ) as discussed in the main text. First notice that putting every 5n(G) to 
zero takes us back to the uniform charge density case considered by Murthy and Shankar 
f2~5fl. They showed that the canonical transformation in that case is 



f / (A)^e^ = e ( A ^^ ct(q)A(q) - h - C -)), (62) 



where 9 = Xfln^Td and A = arctan/z//z with /x 2 = l/lu — 1. The value of the constant A is 
fixed by requiring that the Hamiltonian in the FR doesn't have a term coupling the particle 
and the oscillator degrees of freedom. 

The Hamiltonian (|PT| ) is different from the one considered in [^J by having terms pro- 
portional to 5n/n. The same is true for the commutator of the kinetic momenta flT9|). We 
will assume that 5n/n is a small parameter and when diagonalizing the Hamiltonian we will 
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only keep terms proportional to it. Consistent with this program a reasonable guess for the 
canonical transformation is 

j3 is a constant that has to be determined later by requiring that the first order of the 



coupling term be zero. In (|62"D c(q) is the operator that corresponds to the uniform density, 
while in fl63|) it depends on 5n. Strictly speaking these are two different operators. 

We proceed as in p5[ by determining the operators A(q, A) and c(q, A) in the new 
representation. Each of these operators is a sum of an unperturbed part that coincides 
formally with 5n = result and a first order in 5n/n part. We introduce the notation 
Ao(q, A) + Ai(q, A) for these parts (similarly for c(q, A)). Using the canonical transformation 
(|63|) and the commutation relations (|18|) and (|i~9f ) we derive the following first order flow 
equations for the operators 
dA 1 (q,X) 

"lA Mq,A) 

- P £ c o(q " G > X)HG)UQ^G) +} (64) 

G 

dci(q, A) 



dA 



2eS*nMi(q, A) + 9 (2eB*((3n + 1) - 47r/n) 



x £ A (q - G, A)5n(G)g_(g - (65) 

G 



Substituting Eq. (|65f) into Eq. (p4|) a second order inhomogeneous ordinary differential 
equation is obtained for Ai(q, A). The general solution depends on two arbitrary constants 
that are determined through the initial conditions Ai(q, 0) = ci(q, A) = 0. The result of the 
calculation is 

Ai(q, A) = £ - G)5n(G)g_(g^G) + 

+ (o m) sm ^ A o~Y^~ E c (q - G)5n(G)g_(g - (66) 



'2/in /ieS* ' 2/x 2 



G 



where a new constant a = 2fi 2 (/3/9 + nl/eB* — l/2n) was introduced. Having found Ai(q, A), 
we can integrate Ci(q, A) from the Eq. (|65|) with the result 
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a A sin fi\ 
2 V G 



ci(q,A) = -=^E<q-G)fc(G)i(g^) + 



+ (<& " ^ *»* + ^) f ^ - ^(G)^ - ^ (67) 
The only term in the Hamiltonian (|TT| ) that is not expressed through the operators A(q) 
and c(q) is the CP kinetic energy T = X^n^^m. We will compute this operator in FR 
by deriving first the flow equation for it. First though we have to rearrange T, by using the 
canonical momentum commutator, into 

m ^IL-IL+ "feB* Til . \ 

T = ^^ + ^{-^-^ ) )- (68) 
The second term in fl68|) will not contribute to the flow equation after applying the RPA 



approximation to it. It turns out to describe the magnetic moment of the CP. After doing 
the appropriate commutators we find that to first order in 5n/n the kinetic energy operator 
T\ obeys the flow equation 

c^A) = eB*l £ (A j (q> A)co(q; A) + ^ A)ci(q; A) + h . c .) 

„ r>* a O-rrlf) 5 5 

+— E E W(q, A)c (q - G, \)6n(G)Ul ~ G)+ + h.c.) . (69) 

m q G 

We can integrate the kinetic energy from (|69|) using the initial condition Ti(0) = 0. The 
resulting expression for the kinetic energy in FR is substituted together with the operators 
A(q, A) and c(q, A) in FR into the Hamiltonian (|I7D . We fix the constant (3 by requiring 
that there be no coupling between the particle and the oscillator degrees of freedom. That 

way we get 

[i + (fi 2 - 1) arctan/i 

/3 = -- — y -. (70) 

4n 2 fj, 2 v 7rZ arctan // 

The other consequences of this transformation are given in the main text. 

IX. APPENDIX II 

To calculate the matrix elements (ni \ p (q) \ri2) and (ni \ pi(q, G) |n 2 ) one needs to know 
what the corresponding matrix elements for the operators e* q r , qxll*e ,q ' r and qxll*e* ( ' q_c; - ) ' r 
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are. We will give here only the final expressions for these matrix elements, as the first two 
were derived in several papers (see for example reference ||) and the third can be found 
using a similar approach. The matrix elements are 



e jqr \n 2 ) 



If (ni| q x IIV q - r |n 2 ) 



= ^\\T2 {qx + iqy) V 

I j" / • \ ni-ri2 



+ K + l)^-r(^o 2 /2))e 



-g 2 J * 2 /4 



If (m|qx n*e^ q - G > r \n 2 ) 




n 1 -n 2 



ni! Vv^ 
q (q-G) 



Ll\- n > ((q-G) 2 /S 2 /2) 



V (q-G) 2 

-n x I%-? ((q-G) 2 /* 2 /2 
+(n 2 + l)^;r ((q-G) 2 /* 2 /2)) + 



qx G 
(q-G) 2 



(71) 



(72) 



(ni - n 2 ) 



v r ni— n2 



q-G) 2 /* 2 /2 e 



-(q-G) 2 iS 2 /4 



(73) 



X. APPENDIX III 

In this Appendix we will illustrate the calculation of the integrals over the momentum q 
that appear in the HF Hamiltonian Eq. (|39l)-(f46l). As an example we will take the integral 
that appears in the exchange contribution of the two-body, first order in 5n term. Other 
integrals are done in a similar way. We choose to integrate the following term which is part 
ofEq. 

Bifa=^£V(q) (0|Po(-q)|l> 

x (0| pi(q, G) |0> e -«S 2 qx(G+ 2Gl )/2_ (74) 

Using the formulas given in the Appendix II for the density operator matrix elements we 
find 
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27re 2 g f°° dq 



2 TV 



U tbe = — ^ / -\ 2 / d6e-i 
2e Jo {2tc) 2 Jo 



A 2 



X c 



,q G c 



G 2 2(c+l) 



xe -(q-G) 2 /S 2 /4 e -^S 2 qx(G+2G 1 )/2 > ^ 

Taking into account that q • G = q(G^e ie + G+e- t9 )/2 and q x G = q(G + e- w - G^e i9 )/2i, 
where G + = G x + iG y and G_ = G x — iG y , we get 

d6 



2eZ5 7o io 2tt 



;/) X 



-id 



x ze — = c 



2y/2) \ 2GHI 

x 2 + x{G-e w + G+e- ie )Q 



2(c+ 1) 

xe ((G„+G 1 _)e«-G 1+ e-^)«Sx/2 e -Z S 2 G 2 /4 ) ^ 

where x = ql$. First we will integrate with respect to the variable x. We notice that it is 
possible to extend the interval of the integration over the whole real axis. The integrand 
in Eq. ( |76D is such that the odd/even powers of x are multiplied by the exp(i9n) with n 
odd/even. Then reversing the sign of x and making a transformation 9' = 9 + 7r doesn't 
change the integrand while shifting the integration with respect to x interval to (— 00, 0). 
The integral is then 

ie 2 g d9 ( (G- + G + e~ 2ie ) 



the 



l* Jo 2tt V 4V ,HI 4V2 



c 2 cln , r , ce 10 \ _ z .2 G 2 /4 



x( ^-^ ) - /5( ^ ) 47^J e "^ (77) 

We introduced a notation for the Gaussian integral I n (a,/3) = Jf^ dx exp(—ax 2 + 2f3x)x n 
with a = A 2 /l* 2 + 1/2 and (3 = ((G_ + Gi_)e ie - G 1+ e~ i0S \ l* /A. An important observation 
is that after Eq. ([77]) is expanded the result is the sum of the integrals of the form (note 
the even powers of e %e that appear) 



Jo 2n a 
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2w (ae^ pl2n9 de I ifn>0 

' rrr otherwise, 

|n|! ' 



2 - ru„-ie,2 d# if n < 







,(fee-' 9 ) 2 r ,i2rae uu 
27T 



i.! 



(78) 
otherwise. 



Here a = (G_ + and 6 = — G\ + . The second and third lines in Eg. ([78]) are given 

because they are used to calculate other integrals. They hold if either one of a or b are zero. 
The final answer is then a series of the modified Bessel functions multiplied by appropriate 
constants. 



35 



REFERENCES 

[1] D. C. Tsui, H. L. Stormer and A. C. Gossard, Phys. Rev. Lett. 48, 1599, (1982). 
[2] R. B. Laughlin, Phys. Rev. Lett. 50, 1395, (1983). 

[3] Composite Fermions: A Unified View of the Quantum Hall Regime, edited by 
O. Heinonen, (World Scientific, Singapore, 1998). 

[4] J. K. Jain, Phys. Rev. Lett. 63, 199 (1989); Phys. Rev. B 41, 7653 (1990). 

[5] S. M. Girvin, in Chapter 9 of "The Quantum Hall Effect", Edited by R. E. Prange and 
S. M. Girvin, Springer- Verlag, 1990; S. M. Girvin and A. H. MacDonald, Phys. Rev. 
Lett. 58, 1252 (1987); S. C. Zhang, H. Hansson and S. A. Kivelson, Phys. Rev. Lett. 
62, 82 (1989); N. Read, Phys. Rev. Lett. 62, 86 (1989); D. H. Lee and S. C. Zhang, 
Phys. Rev. Lett. 66, 1220 (1991). 

[6] A. Lopez and E. Fradkin, Phys. Rev. B 44, 5246 (1991); ibid 47,7080 (1993); Phys. 
Rev. Lett. 69, 2126 (1992); V. Kalmeyer and S. C. Zhang, Phys. Rev. B 46, 9889 
(1992); B. I. Halperin, P. A. Lee and N. Read, Phys. Rev. B 47, 7312 (1993). 

[7] R. Shankar and G. Murthy, Phys. Rev. Lett. 79, 4437 (1997). 

[8] G. Murthy and R. Shankar, Phys. Rev. B 59, 12260 (1999); G. Murthy, Phys. Rev. 
B 60, 13702 (1999); R. Shankar, Phys. Rev. Lett. 83 2382 (1999); G. Murthy, Jour. 
Phys. Cond. Mat. 12, 10543 (2000); R. Shankar, Phys. Rev. B 63 085322 (2001). 

[9] See articles by H. A. Fertig and M. Shayegan in Perspectives in Quantum Hall Effects, 
edited by S. Das Sarma and A. Pinczuk, (Wiley, New York, 1997). 

[10] K. Maki and X. Zotos, Phys. Rev. B 28, 4349 (1983). 

[11] D. Levesque, J. J. Weiss and A. M. MacDonald, Phys. Rev. B 30, 1056 (1984). 
[12] P. K. Lam and S. M. Girvin, Phys. Rev. B 30, 473 (1984). 



36 



[13] H. W. Jiang, R. L. Willett, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer and K. W. West, 
Phys. Rev. Lett. 65, 633 (1990). 

[14] H. W. Jiang, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer and K. W. West, Phys Rev. B 
44, 8107 (1991). 

[15] H. Fukuyama, P. M. Platzman and P. W. Anderson, Phys. Rev. B 19, 5211(1979). 

[16] D. Yoshioka and H. Fukuyama, J. Phys. Soc. Jpn. 47, 394(1979). 

[17] V. J. Goldman, M. Shayegan and D. C. Tsui, Phys. Rev. Lett. 61, 881 (1988). 

[18] E. M. Goldys, S. A. Brown, R. B. Dunford, A. G. Davies, R. Newbury, R. G. Clark, 
P. E. Simmonds, J. J. Harris and C. T. Foxon, Phys. Rev. B 46, 7957 (1992). 

[19] V. J. Goldman, J. K. Wang, B. Su and M. Shayegan, Phys. Rev. Lett. 70, 647 (1993). 

[20] H. Yi and H. A. Fertig, Phys. Rev. B 58, 4019 (1998). 

[21] L. Zheng and H. A. Fertig, Phys. Rev. B 50, 4984 (1994). 

[22] L. Zheng and H. A. Fertig, Phys. Rev. Lett. 73, 878 (1994). 

[23] L. W. Engel, T. Sajoto, Y. P. Li, D. C. Tsui and M. Shayegan, Surf. Sci. 263, 44 (1992). 

[24] H. Fukuyama and P. A. Lee, Phys. Rev. B 17, 535 (1978); B. G. A. Normand, P. B. Lit- 
tlewood and A. J. Millis, Phys. Rev. B 46, 3920(1992); H. A. Fertig, Phys. Rev. B 59, 
2120 (1999). 

[25] G. Murthy and R. Shankar in reference 1. 

[26] Note that quantities of the form V± mean V x ± iV y . 



[27] Equation ( |2"2|) differs from its uniform density counterpart [^SJ because the canonical 
momentum operator that enters the definition of c(q) here depends on the density 
modulations through the CS gauge field. 



37 



[28] W. Kohn, Phys. Rev. 123, 1242 (1961). 

[29] S. M. Girvin. P. M. Platzmann, and A. H. MacDonald, Phys. Rev. B 33, 2481 (1986). 

[30] The operator c(q) depends implicitly on 8n, because IT, contains the vector potential 
A*(r J ), which in its turn is related to the average density through V x A*(i\,-) = B — 
2itln/e — 2nl8n(rj)/e. 

[31] R. Cote and A. M. MacDonald, Phys. Rev. B 44, 8759 (1991). 

[32] D. Yoshioka and P. A. Lee, Phys. Rev. B 27, 4986 (1983). 

[33] When 1 there are overlapping bands so that the bands that overlap with Mg'th 

and Mq + l'st bands should be included in determining the excitation gap. 

[34] D. S. Fisher, B. I. Halperin and R. Morf, Phys. Rev. B 20, 4692 (1979). 

[35] E. Cockayne and V. Elser, Phys. Rev. B 43, 623 (1991). 

[36] We looked into the possibility that the quasihole lattice may be formed but found that 
its ground state energy is higher in the relevant range of filling factors. 

[37] This phenomenon is analogous to the peak effect in superconductors. See, for example, 
G. Blatter, M. V. Feigelman, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Rev. 
Mod. Phys. 66 1125 (1994). 



38 



- 


□ 


- 




□ 

* * 






□ 

* □ 






□ 








* * 

* 




> 


□ * 



0.16 0.18 0.2 0.22 0.24 

V 

Figure 1. 

Caption: The activation gap dependence on the filling factor around v « 1/5. Squares 
are our CF theory with four vortices attached for the hexagonal lattice. Diamonds represent 
our CF theory with four vortices attached for the oblique lattice (see the text). Stars are 
experimental results read off Fig. 3 of reference ||14j| . 
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Figure 2. 

Caption: Effective HF potential for different values of parameter A. The filling fraction 
specific factor (1— lv) 2 was omitted from the zeroth-order expression of the effective potential 
when generating these curves. 
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Figure 3. 

Caption: The activation gap dependence on the filling factor around v fa 1/5 (/ = 2). 
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Figure 4. 

Caption: The shear modulus fi for the triangular CF lattices as a function of filling factor 
(i/ < 1/5). 
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Figure 5. 

Caption: The HF energies for CF lattices differing by an angle 9 between the reciprocal 
lattice vectors for the filling factor v ~ 0.206. The zero on the vertical axis corresponds to 
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Figure 6. 

Caption: The activation energies for CF lattices differing by an angle 9 between the 
reciprocal lattice vectors for the filling factor v ps 0.206. 
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